Spin Hall effect in graphene due to random Rashba field 
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Spin Hall effect due to random Rashba spin-orbit coupling in the two-dimensional honeycomb 
lattice of carbon atoms (graphen) is considered theoretically. Using the Green function method and 
diagrammatic technique we show that fluctuations of the Rashba interaction around zero average 
value give rise to nonzero spin Hall conductivity. Generally, the conductivity is not universal, but 
depends on the ratio of the total momentum and spin-flip relaxation rates. 
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Introduction: Interest in the phenomena due to spin- 
orbit coupling has revived in recent years owing to the 
possibihty of pure electrical control of the spin degree 
of freedom. One of such phenomena is the spin Hall ef- 
fect (SHE), predicted by Dyakonov and Perel in 1971—"—. 
The effect consists in a spin current (or spin accumula- 
tion) flowing perpendicularly to an external electric field 
in systems with strong spin-orbit interaction. Of partic- 
ular interest is the SHE in nonmagnetic systems, where 
the induced spin current is not accompanied by a charge 
current. Physical mechanisms that lead to SHE arc ei- 
ther of intrinsic or extrinsic origin (for review see refs 
0-0) ■ The intrinsic SHE is associated with band struc- 
ture modified by an intrinsic spin-orbit coupling of the 
material under consideration, and is a consequence of an 
unusual trajectory of electrons in the momentum space. 
Extrinsic contribution to the SHE, in turn, is a conse- 
quence of spin-orbit interaction due to impurities and 
other structural defects that lead to spin-dependent scat- 
tering - skew scattering and side jump. 

In a general case, spin-orbit coupling in a system con- 
tains regular (spatially uniform or periodic) and random 
components^. Local imperfections of the system, like 
random distribution of donors or impurities, may lead to 
local enhancement/suppression of the regular spin-orbit 
coupling and as a consequence charge carriers in the sys- 
tem propagate in an effective fluctuating spin-orbit field. 
Such fluctuations of spin-orbit field strongly modify spin 
relaxation and significantly affect spin transport. Dugaev 
et al^ reported recently that the random Rashba field 
in two-dimensional electron gas generates nonzero SHE, 
even in the presence of potential scattering by impuri- 
ties. This is distinct from the well known result for a 
two-dimensional electron gas with constant Rashba cou- 
pling, where potential scattering by impurities cancels 
out the universal intrinsic spin Hall conductivityiSi^. 

In this paper we consider SHE created by fluctuat- 
ing Rashba spin-orbit coupling in graphene - a two- 
dimensional honeycomb lattice of carbon atom o^^'^'^ . The 
corresponding intrinsic spin-orbit coupling in graphene is 
very weak^^, so the quantized value of spin Hall conduc- 
tivityi^, that appears inside the spin-orbit energy gap, 
is extremely difficult to be verified experimentally. In 



turn, the Rashba type spin-orbit interaction in graphene 
appears generally due to a substrate and can be rela- 
tively strong. This interaction includes usually a regular 
component, but a random term may also appear due to 
ripples of the graphene plane, disorder, electron-phonon 
coupling in the substrate, presence of adatoms, cto^^—. 
The influence of fluctuating Rashba fleld on spin relax- 
ation in graphene has been analyzed for two different 
models of the Rashba fleld fluctuation a^^'^° . Here we 
show that the spin Hall conductivity associated with such 
fluctuations is not universal and depends on the ratio 
of total momentum and spin relaxation times. We con- 
sider two different types of the correlation function for 
the spin-orbit coupling strength. Below, we briefly de- 
scribe the model under consideration and the theoretical 
method used to calculate spin Hall conductivity. Then, 
we calculate the Feynman diagrams representing various 
contributions to the conductivity and also evaluate the 
spin dephasing time. 

Model and method: Since the intrinsic spin-orbit inter- 
action in a free standing graphene is very small, it will be 
neglected in the following. Accordingly, the low-energy 
electronic states around the Dirac point K are described 
by the Hamiltoniani^ , which in the sublattice space 
takes the matrix form 



K 





v{kx, + iky) 



v{kx - iky) 




(1) 



where A:2-(y) ^^re the wavevector components, and v 
with vp being the electron Fermi velocity. 

As already mentioned above, an important spin-orbit 
interaction in graphene is the Rashba coupling due to a 
substrate - though a related interaction may also appear 
in a free standing graphene with ripples. In the following 
we consider the case when the Rashba coupling vanishes 
on average, but the corresponding coupling parameter 
fluctuates in space around zero. These fluctuations of 
the Rashba field will be treated as a perturbation, so 
the full Hamiltonian for the point K can be written as 
= Hq + V^^,, where the perturbation V^, takes 



the form 



Akk' 
Akk'(0'y - iCTx) 







(2) 
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with Ak k' describing the corresponding Rashba couphng 
term, and cTq, (a = x, y, z) being the Pauli matrices in 
spin space. According to our assumptions, (Aq) = and 
(A^) = Cq^Oforq = k-k'. 

Following the Kubo-Streda formula^ 
spin Hall conductivity in the form 



we write the 



(3) 

The term crj.|,^^ is the so-called topological contribution 
to the spin Hall conductivity, which comes from states be- 
low the Fermi level. This contribution, however, is equal 
to zero for the model under consideration {cr^y^^ ^ only 
inside the energy gap induced by intrinsic spin-orbit in- 
teraction, which has been neglected here). In turn, the 
term a^^^ is determined by states at the Fermi level and 
is given by the formula 

„ , eh f (fk 
2^ 
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i2ny 



Tr{,f,'G"{eF)vyG''{sF)}, (4) 



where G'^^"^'(eF) is the retarded (advanced) Green func- 
tion at the Fermi level e^; jx' = (1/2) [vx, Sz]^ is the spin 
current density operator, and Vi = {l/fi){dH/dki) is the 
velocity operator. In the model under consideration, j^' 
and Vy take the following explicit forms: 



Jx 



CT. 




and 











(5) 



(6) 



where ctq is a unit matrix in the spin space. 

Spin Hall conductivity and spin relaxation: Perturba- 
tion expansion for the Green functions in Eq. (4) leads to 
a series of Feynman diagrams contributing to the spin 
Hall conductivity (see Fig.l). Thus, one can write 
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The term (Jxy represents a contribution from the bare 
bubble diagram and is generally a part of intrinsic spin 
Hall conductivity. Since the nonperturbative part of the 
Hamiltonian does not contain any term associated with 



spin-orbit interaction, this term vanishes, axy 
The second term in Eq.(7) can be written as 



(1) 



0. 



27r 



'xy 



(8) 



(27r)2 (27r)2 
X Tr [j^ Vkk' G^, VyG^, Vk'kGj^] 

In turn, the last term of Eq.(7), corresponding to the 
bottom diagrams in Fig.l, takes the form 



xy 



d^k d^k' 



2n J (27r)2 (27r)2 
xTr [.f/G^VyG^Vuu'G^,Vu'i.G^ 

+J^^G^Vkk'G^.Vk'kG^,Gl^] . (9) 





FIG. 1. Feynman diagrams corresponding to various contri- 
butions to the spin Hall conductivity in the dc limit. 



The Green functions G^'^^^ have the form, G^^^^ = 

[{ep ± iT) - H^] ' \ where F = h/2T. We assume that 
T is the total relaxation time, which includes scattering 
from impurities, r^, and scattering from random Rashba 
field, Tg. Thus one can write 1/t = l/ri + l/Ts- Note that 
for small impurity density the scattering from random 
Rashba field is the dominant mechanism of relaxation 
{t ^ Tg). To obtain we calculate the self energy, 

|J^lik'G°«Fk,k. (10) 

Taking into account the explicit form of the Green func- 
tion G'^^ corresponding to the Hamiltonian (1), one ob- 
tains 



d^k' 



(27r)2 (ep ^ Elk' +iO+)ieF - E2k' 
where i?i 2 = ±vk = isk, and 



M=( ^('^o -cT.) 
U 9(0-0 



-iO+)' 
(11) 



(12) 



^From the above equations follows that the self energy 
can be written as 



Zn 



de Cq 5{eF - fk-q) 



-iM- 



2kf 



dq 



2Ca 



or equivalently 



(13) 



where F^ = Tl/2ts- From Eq.(13) follows that the explicit 
form of Tg depends on the form of the correlator G^iS^. 
Let us go now back to Eq.(8), which can be rewritten 

as 

e r d'^k d'^k' 
2^ J (2^(2^ 
l6\Xkk'\^v\kxk'^-kyk'y)eFr 
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{eF + iTY){el,-{eF+iTYy 



(14) 
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s (2) 

Evaluation of this formula gives axy = 0. Thus, the 
spin Hall conductivity is given by the bottom diagrams 

in Fig.l, a'^'y = cr^y^'^^'*\ and for weak disorder one finds 



a. 



e V Sp 



xy 



dq q / dk k 
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X / deC,5{sl~e%)5isl_^~s%). (15) 



Finally the spin Hall conductivity can be written in the 
following simple form: 



2k f 



Co 



V^^/Ikj^~(f 



(16) 



Taking into account the formula for spin relaxation 
time, and including contributions from both K-points, 
we may write the formula for spin Hall conductivity in 
the form, 



'xy 



e T 
47r r. 



(17) 



This formula clearly shows that SHE is generally not 
universal. Potential scattering by impurities reduces the 
spin Hall conductivity, but does not suppress it to zero. 

The spin Hall conductivity depends on the form of the 
correlator Cq. Dugaev et al^^ assumed the correlator in 
the form 



^(1) 



(18) 



Such fluctuations may originate from the ripples or some 
impurities that exist at the surface of graphene. Using 
this correlator one obtains (for a single K point) 
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FIG. 2. Spin Hall conductivity as a function of kpR for 
constant (as indicated) values of R. The parts a) and b) cor- 
respond to the correlation function defined by Eq.(17) and 
Eq.(20), respectively. Contributions from both K points are 
included. The parameters used in the calculations are as 



follows: (A^) = 25 X lO"-* eV^ v = 3.516 x 10" 



eVm, 



n = 3 X 10' 



and r = 6.58 x 10"** eV. 



<l = 7f^(A')fcFi?' {WkpR) - Lo{2kFR)) , (19) 

where /q and Lq are the modified Bcsscl and Struve func- 
tions of zeroth order, respectively. From this relation one 
finds asimple formula in some limiting cases. 



47r vV 



R 



TrkpR ; kpR < 1 
1 ;fcFi?>l 



(20) 



In turn, Zhang and Wui^ have proposed recently an- 
other form of the correlator that may describe fluctua- 
tions of the Rashba field due to the presence of adatoms 
distributed on top of the graphene surface and between 
the graphene layer and substrate. 
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4:TT^n{X^)R^e 



(21) 



where n is the impurity density. The corresponding spin 
Hall conductivity (for single K point) takes the form 



and in the limiting cases 



{\^)kFR'^nIo{2kj.R^)e 
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en 
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nkpR -^kpR < 1 
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■^kpR > 1 



(22) 



(23) 



Numerical results: Numerical results on the spin Hall 
conductivity in a general situation (arbitrary kpR) are 
presented in Fig. 2 and Fig. 3. The contributions from 
both K and K' points are included there. Figure 2 shows 
the spin Hall conductivity as a function of kpR for fixed 
values of the correlation length i?, as indicated. Parts 
(a) and (b) correspond to the two considered forms of 
the correlator Cq. In turn, Fig. 3 shows the spin Hall 
conductivity as a function of kpR in the case when kp is 
fixed as indicated. As before, the upper and lower parts 
correspond to the two forms of the correlator Cq. 

From these two figures clearly follows that the spin Hall 
conductivity vanishes when either the Fermi wavevector 
kp or the correlation length R tend to zero. To un- 
derstand this behavior one should bear in mind that 
kp = corresponds to the Fermi level at the Dirac 
points, where the density of states vanishes. In turn, the 
limit of i? = corresponds to the electron Fermi wave- 
length much longer than the scale length of the disor- 
der. These limiting values of the spin Hall copnductivity 
are in agreement with the corresponding formulas (20) 
and (23). When kpR grows (keeping R constant), the 
conductivity initially grows with kpR and after reaching 
maximum it tends to the limits for kpR ^ 1 according 
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FIG. 3. Spin Hall conductivity as a function of kpR for 
constant (as indicated) values of kp- The parts a) and b) 
correspond to the correlation function defined by Eq.(17) and 
Eq.(20), respectively. Contributions from both K points are 
included. The other parameters are as in Fig. 2. 



to the formulas (20) and (23). Note, these limits depend 
on R, as follows from (20) and (23). In turn, when kp 



is fixed (see Fig. 3), the spin Hall conductivity increases 
with increasing kpR, but the rate of this increase depends 
on fci^. 

To conclude, we have analyzed the spin Hall effect in- 
duced by fluctuating Rashba spin-orbit interaction in a 
monolayer graphene. The key assumption was the ab- 
sence of uniform Rashba term, i.e. the assumption that 
the Rashba interaction vanishes on average. Two dif- 
ferent forms of the correlation function have been con- 
sidered, which may correspond to different physical sit- 
uations. We note, that such random Rashba interac- 
tion may appear as a result of ripples in a free stand- 
ing graphene or due to impurities on both sides of the 
graphene plane. 

It has been shown, that the spin Hall conductiv- 
ity induced by the fluctuating Rashba held is equal to 
(e/47r) (t/ts), independently on the form of the correla- 
tor, and is generally not universal. A universal value oc- 
curs only when both r and are equal (as in the absence 
of any other defects, except the fluctuating Rashba term) . 
It is thus evident that potential scattering by defects sup- 
presses the spin Hall effect, but this suppression is not 
complete, similarly as in the case of two-dimensional elec- 
tron gas. 
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